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We study the quantum communication protocol of remote state preparation (RSP) for pure states 
of qubits encoded in single photons transmitted through a double slit, the so-called spatial qubits. 
Two measurement strategies that one can adopt to remotely prepare the states are discussed. The 
first strategy is the well-known spatial postselection, where a single-pixel detector measures the 
transverse position of the photon between the focal and the image plane of a lens. The second 
strategy, proposed by ourselves, is a generalized measurement divided into two steps: the imple- 
mentation of a two-outcome positive operator- valued measurement (POVM) followed by the spatial 
postselection at the focal plane of the lens by a two-pixel detector in each output of the POVM. In 
both cases we analyze the effects of the finite spatial resolution of the detectors over three figures 
of merit of the protocol, namely, the probability of preparation, the fidelity and purity of the re- 
motely prepared states. It is shown that our strategy improves these figures compared with spatial 
postselection, at the expense of increasing the classical communication cost as well as the required 
experimental resources. In addition, we present a modified version of our strategy for RSP of spatial 
qudits which is able to prepare arbitrary pure states, unlike spatial postselection alone. We expect 
that our study may also be extended for RSP of the angular spectrum of a single-photon field as an 
alternative for quantum teleportation which requires very inefficient nonlinear interactions. 

PACS numbers: 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

Remote state preparation (RSP) [1-3] is a quantum 
communication protocol similar to quantum teleporta- 
tion (QT) [4], regarding their ultimate goals: a sender, 
Alice, wants to transmit to a receiver. Bob, the quan- 
tum state of a qubit without sending physically the in- 
formation carrier. In both protocols Alice and Bob need 
to share entanglement and the state of Bob's system is 
remotely prepared, after Alice has communicated classi- 
cally to him the outcome of a measurement on her sys- 
tem. However, unlike QT, in RSP Alice has full knowl- 
edge of the state to be sent. For this reason, while in 
QT the consumption of one maximally entangled state 
(ebit) and two bits of classical communication (cbits) are 
both necessary and sufficient for faithfully transmitting 
the state, in RSP it is possible to trade off between these 
two resources [1-3, 5, 6]. For instance, a qubit lying in 
a given great circle of the Bloch sphere can be remotely 
prepared with 1 cbit from Alice to Bob, if they share 1 
ebit of entanglement [2]. Also, it is possible that Alice 
and Bob share a nonmaximally entangled state (i.e, less 
than 1 ebit) and she can remotely prepare his qubit ei- 
ther in a state within specific ensembles using 1 cbit or 
arbitrary pure states using finite (more than 1) cbits [6]. 

Another important aspect that distinguishes RSP from 
QT is the measurement strategy which Alice has to em- 
ploy to accomplish the protocol. In QT Alice needs to 
perform a complete Bell-state measurement (BSM) be- 
tween two particles, one from the entangled pair shared 
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with Bob and one carrying the qubit encoding the state 
to be transferred. For optical implementations, this is a 
problem because complete BSM is not possible with lin- 
ear optics only [7]. In RSP, however, the measurement 
strategy is more flexible. In the simplest case it may only 
require projective measurements, as for instance, when 
Bob knows in advance that Alice will only prepare states 
lying in a great circle of the Bloch sphere [2] . When the 
resources for RSP are the same as for QT (1 ebit, 2 cbits), 
Alice will have to implement a generalized measurement 
on her qubit. This kind of measurement, which is treated 
within the formalism of positive operator-valued mea- 
surement (POVM), can be implemented through a uni- 
tary interaction between the Alice's system and an an- 
cillary one, followed by measurements of both systems 
[8]. The ancillary system may be just a second degree of 
freedom of Alice's photon and thereby there is no con- 
straint for implementing the POVM with linear optics 
only [9, 10]. 

So far, several optical realizations of RSP of single 
qubits have been reported, including pure state of qubits 
in a superposition of single-photon and vacuum states 
[11], pure and mixed states of photon polarization qubits 
[12-16], and pure states of spatially encoded photonic 
qubits [17, 18]. These realizations were either probabilis- 
tic [11-13, 17, 18] or deterministic [14-16]. In the first 
case, which some authors refer to RSP just as conditional 
preparation [19], there is a cost of 1 cbit and the unitary 
correction by Bob is not required. In the deterministic 
realizations, [14, 15] have successfully implemented the 
POVM but not the required unitary correction and there- 
fore constituted a proof-of-principle demonstration of the 
RSP protocol. In Ref. [16], however, the authors have 
successfully implemented both the POVM and Bob's uni- 
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tary correction, making the full demonstration of RSP 
with a cost ranging from cbits for preparation of max- 
imally mixed state to 2 cbits for arbitrary pure states. 

In this article we study the remote preparation of pure 
qubits states encoded in the discretized one-dimensional 
transverse momentum of single photons. These so- 
called spatial qubits are obtained when a paraxial 
and monochromatic single-photon field is made to pass 
through a double slit aperture [20]. The advantage of- 
fered by this kind of encoding is that one can create 
high-dimensional qudits simply by increasing the num- 
ber of slits in the aperture [20, 21]. Recent experimental 
works, employing photon pairs from spontaneous para- 
metric down-conversion, have demonstrated the entan- 
glement between two spatial qubits [17, 22-24] as well 
as for two spatial qudits [21, 25]. Among other applica- 
tions, spatial qubits and qudits have been used to demon- 
strate hybrid photonic entanglement [26] and a variety of 
quantum tomographic techniques as single scan tomog- 
raphy [17, 25] and tomography based on measurements 
either onto mutually unbiased bases [27, 28] or symmet- 
ric POVMs [29]. Those applications have benefited of 
the considerable effort for finding methods to measure 
spatial qubits and qudits [17, 21, 22, 25, 28, 30, 31]. 

Here, we assume that Alice and Bob share 1 ebit and 
we discuss two measurement strategies that she could 
adopt to remotely prepare his spatial qubit. The first 
strategy is the well-known spatial postselection, where 
a single-pixel detector measures the transverse position 
of the photon between the focal and the image plane 
of a lens [17, 21, 22, 25]. The second strategy, some- 
what similar to the one proposed by ourselves in [30], 
is a generalized measurement divided into two steps: (i) 
the implementation of a two-outcome POVM where the 
photon polarization is used as the ancillary system; (ii) 
the spatial postselection at the focal plane of the lens 
by a two-pixel detector in each output of the POVM. 
Both strategies share the common feature that the ul- 
timate physical property to be measured is the photon 
transverse position. In this case a realistic description 
demands us to take into account the finite spatial res- 
olution of the detection system, where the best it can 
be done is to locate the photon within a narrow interval 
rather than a point. For both, spatial postselection and 
generalized measurement, we analyze, numerically, the 
effects of a realistic detector over three figures of merit 
of the RSP protocol, namely, the probability of prepara- 
tion, the fidelity and the purity of the remotely prepared 
states [32]. It is shown that our strategy improves these 
figures compared with spatial postselection, at the ex- 
pense of increasing the classical communication cost as 
well as the required experimental resources. In addition, 
we briefly discuss the RSP of spatial qudits and show 
that by using the generalized measurement strategy we 
proposed in [30] one can prepare arbitrary pure states, 
unlike spatial postselection alone. 

Before proceeding, it is worth to draw a comparison 
between RSP of spatial and polarization qubits. In the 



latter, a probabilistic RSP of arbitrary pure states is re- 
alized with 50% efficiency through projective measure- 
ments [13]. On the other hand, a deterministic imple- 
mentation employs POVMs and it is 100% efficient (in 
principle) [14-16]. While it is obvious that the prob- 
ability of preparation improves, the fidelity and purity 
of the remotely prepared polarization qubits are limited, 
in both cases, only by the precision of the measurement 
instruments (polarizers, wave plates, etc.) so that they 
could be one in principle. For spatial qubits, the strategy 
of spatial postselection would be the analog of projec- 
tive measurement in polarization. However, here the ef- 
ficiency is smaller than 50% since the photon propagates 
in a continuous space confined within the diffraction en- 
velope, and the detector measures only a small fraction 
of this space. Similarly, our proposed generalized mea- 
surement strategy would be the analog of a deterministic 
implementation for polarization qubits, but, again, due 
to the diffraction we are not able to get 100% efficiency. 
Therefore, unlike polarization, for spatial qubits the re- 
lation between the figures of merit is not trivially estab- 
lished when one compare the two strategies, since these 
figures are affected by the diffraction and/or the finite 
spatial resolution of the detector. Those are the aspects 
we will approach here. 



This article is organized as follows: in Sec. II we out- 
line the general framework for RSP of spatial qubits. 
In Sec. Ill we discuss the spatial postselection strategy 
starting with the idealized description of point detectors 
(Sec. HI A) and then moving on to the realistic descrip- 
tion of detectors with finite spatial resolution (Sec. HI B) . 
In the latter case we derive expressions for the three fig- 
ures of merit to be analyzed and compare our theory with 
the experimental results from [17]. A numerical analysis 
of the figures of merit is presented in Sec. HI C. Section IV 
describes the generalized measurement strategy proposed 
in this work, first presenting the POVM (Sec. IV A) and 
then the detection system. Here, we also begin by consid- 
ering the idealized point detectors (Sec. IV B) and next, 
describe the realistic case of detectors with finite resolu- 
tion and their effects on RSP (Sec. IV C). We evaluate 
the figures of merit in this latter case and present our 
numerical results (Sec. IV D). In Sec. V we compare the 
two strategies regarding the figures of merit as well as 
the resources required to implement the protocol. Sec- 
tion VI describes the RSP of spatial qudits either through 
spatial postselection or the generalized measurement pro- 
posed in [30] , showing the advantages of the latter over 
the former. Finally, in Sec. VII we conclude the arti- 
cle and briefly discuss a possible extension of our results 
for the remote preparation of the angular spectrum of a 
single-photon field as an alternative route for quantum 
teleportation. 
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FIG. 1: (Color online). Schematic representation of RSP of 
spatial qubits. A SPDC source generates spatially correlated 
photon pairs whose transverse momentum is discretized by 
double slits. The two-qubit maximally entangled state (1) is 
then distributed between Alice and Bob. To remotely prepare 
the target state (2), she measures her photon employing one 
of the two strategies to be described here and communicates, 
through a classical channel, her outcome to Bob. Having this 
information, he may or may not perform a unitary correction, 
U , and achieve the desired target state (2). 



II. RSP OF SPATIAL QUBITS 

Entangled states of spatial qubits are created by letting 
each member of the photon pair generated by sponta- 
neous parametric down-conversion (SPDC) cross a dou- 
ble slit aperture as sketched in Fig. 1. The degree of 
entanglement between them can be controlled by manip- 
ulating the pump beam transverse profile [20, 22, 23]. 
Let us assume that this source creates the maximally en- 
tangled state 



(a) 
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(1) 



which is shared between Alice {A) and Bob {B). Here, \r) 
denotes the state of the photon transmitted through 
the right (left) slit. Experimentally, this can be achieved 
by focusing the pump profile at the double slit plane [20, 
22]. Now, Alice's goal is to remotely prepare Bob's qubit 
in the pure state 



\4^)B^a\l)+P\r), 



(2) 



where |a|^-|-|/3|^ = 1. For doing so, she measures her spa- 
tial qubit through one of the two methods to be discussed 
here and communicates to Bob, through a classical chan- 
nel, the outcome of her measurement (see Fig. 1). With 
this information. Bob may or may not have to perform 
a unitary transformation iU) to recover the target state 
given by Eq. (2). This, as we will see, depends on the 
measurement strategy adopted by Alice and the outcome 
of this measurement. 



III. RSP OF SPATIAL QUBITS THROUGH 
SPATIAL POSTSELECTION 

Firstly, let us consider the RSP protocol when Alice 
remotely prepares the spatial qubit in Bob's laboratory 
through spatial postselection. This method relies on the 
detection of the photon transverse position after it has 
propagated through a given optical system. Typically, 
this optical system is comprised of a thin convergent lens 
placed at some distance from the double slit, followed 
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FIG. 2: (Color online), (a) Schematic representation of the 
measurement strategy called spatial postselection. A lens of 
focal length / is placed at a distance 2/ from the double slit. 
A single pixel detector is allowed to move between the focal 
{z = /) and the image {z — 2/) planes in the x and z direc- 
tions. For a point detector, the corresponding measurement 
operator is given by Eq. (6) . For a detector with finite resolu- 
tion (as shown in the inset) the measurement operator is given 
by (14). The regions I and II set the validity of the Fraun- 
hofer and Fresnel approximations, respectively. The former is 
valid if the parameter ka^ jlZ is close to zero as shown by the 
graphic in (b). In the region I the measurement operators are 
calculated with the amplitudes given by (5) while in region II 
they are calculated from (4). 



by a detector which can move in the transverse {x) and 
longitudinal (z) directions between the focal plane of this 
lens and the image plane of the double slit [17, 22, 23, 31]. 
Figure 2(a) outlines Alice's apparatus to implement spa- 
tial postselection. In the following, we discuss two cases 
regarding Alice's detector system: first, the idealized case 
of a point detector and second, the realistic case of a sin- 
gle pixel detector with finite resolution. 



A. Idealized measurement: The point detector 

1. Measuring a spatial qubit with a point detector 

A point detector is, obviously, an idealized device 
which delivers infinite spatial resolution in the measure- 
ment of the photon transverse position. Nevertheless, it 
is important to understand how the RSP protocol would 
work if Alice had such a device at her disposal, before we 
move on to the realistic case of detectors with finite reso- 
lution in the next subsection. Here, we describe the effect 
of a measurement by a point detector over an arbitrary 
state, p, of a spatial qubit. 

Figure 2(a) shows that after transmission through the 
double slit, which discretizes its state, the photon propa- 
gates in the continuous space through the optical system. 
A lens of focal length / is placed at a distance 2/ from 
the double slit. Therefore, from our previous discussion, 
the point detector will be constrained to move, in the 
transverse direction x, between z — f (the focal plane) 
and z = 2f (the plane where the double slit image is 
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formed with a magnification of one). A measurement 
in a given position (x, z) postselects the photon in the 
(non-normalized) state [17] 



\if{x, z)) = Lpi{x, z)\l) + ipr{x, z)\r), 



(3) 



where, up to irrelevant phase factors, the probability am- 
plitudes (pj{x, z) for i = l,r are given by 



pa 



ipj{x,z) 



In this expression a is the slit half-width, 5r — d/2 and 
5i = —d/2, where d is the center-to-center separation be- 
tween the slits and k is the photon wave number. The 
parameters Z — [2p — fz)/{z — /) and rj = {z — f)/f, 
represent the effective longitudinal distance that the pho- 
ton propagates from the double slit to the detection plane 
and a scale factor in the detector's transverse position, re- 
spectively. Equation (4) is the Fresnel diffraction integral 
for a single slit [33] and it is valid everywhere within the 
range f ^ z < 2f where the detector moves. However, 
it will have analytical solution only if the approximation 
Z ^ ka? /2 is satisfied. In this case the quadratic phase 
fmiction in the integral satisfies e^^^ « 1 and Eq. (4) 
reduces to 



ipj{x,z) 
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(5) 



where sinc(^) = sin(^)/^. This is the well known Fraun- 
hofer approximation and Eq. (5) yields the Fraunhofer 
diffraction pattern for a single slit [33]. 

Figure 2(b) shows how the parameter ka^ /2Z changes 
between z = f and z = 2f. This parameter is the maxi- 
mum value of the argument in the quadratic phase func- 
tion in Eq. (4) and determines the regions where the 
Fraunhofer approximation is valid or not. To plot this 
curve we used the values of the parameters a, A and / 
shown in Table I. When the approximation is valid we 
are in the region of Fraunhofer diffraction (region "I" 
in Fig. 2) and the probability amplitudes ipj{x,z) are 
calculated from Eq. (5). Otherwise, we are in the re- 
gion of Fresnel diffraction (region "II" in Fig. 2) and the 
probability amplitudes are calculated from Eq. (4). It 
should be stressed that in both expressions, there is a 
constant term for each plane z which defines the scale 
of the diffraction pattern. This term, together with the 
diffraction envelope itself, has important consequences 
on the efficiency with which Alice can remotely prepare 
Bob's spatial qubit. This point will become clear later. 

From the discussion above the measurement operator 
associated with the detection by a point detector at the 
position (x, z) is [34] 



TABLE L Parameters used in the simulations presented in 
this work: a and d denote the slits half width and the center- 
to-center separation between them, respectively, A is the pho- 
ton wavelength, / is the focal length of the lens, and Ax is 
the detector width. 
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where \ip{x,z)) is given by Eq. (3). The corresponding 
detection probability density for a spatial qubit in a state 
p will be given by (under the assumption that the detec- 
tor has unit quantum efficiency) 



q{x, z) — Tr [s{x, z)p\ 



(7) 



Identifying each point (x, z) with a normalized vector 
of the spatial qubit, Taguchi et al. [17] showed that one 
can access the whole Bloch sphere of the spatial qubit. 
Figure 3 shows some samples of the state vectors one can 
measure by scanning the point detector in x for distinct 
longitudinal planes z. To calculate the vectors for each 
point (x, z) we replaced the probability amplitudes given 
by Eq. (5) into Eq. (3) and normalize them in order to 
obtain their polar and azimuthal angles {9 and 0, respec- 
tively) of the Bloch representation. In addition, in the 
focal plane we considered x G [~T:f/kd,TTf/kd] which 
suffices for a complete loop around the equator, while for 
the other planes z, x took values which ensure that the 
trajectories traced out by the vectors were going from one 
to another pole (see Fig. 3) [35]. The values of the pa- 
rameters a, d, A and / used in this calculation are shown 
in Table I. As shown by the red circles in Fig. 3, at the 
focal plane (z = /) one can measure any state lying only 
in the equator of the Bloch sphere [22, 31]. On the other 
hand, at the image plane (z = 2/) one measure only the 
states {10 7 1^)} (not shown in Fig. 3) which lie at the poles 
of the sphere, what is quite intuitive. The interesting as- 
pect noted by Taguchi et al. [17] is that at intermediate 
planes {f < z < 2f) one can measure arbitrary superpo- 
sitions between \l) and |r). Besides, it is clear from Fig. 3 
that there are several vectors which can be measured at 
more than one plane. This happens when the trajectories 
traced out by the vectors for distinct detection planes z 
intersect each other. For instance, with the exception of 
the focal plane, the logic basis may be measured at any 
plane z € (/, 2/] . In the same way, with the exception 
of the image plane, the state \l) + \r) can be measured 
in X — for any plane z G [/, 2/). Despite this, there 
are differences concerning the efficiency with which these 
"repeated" vectors can be measured, due to the diffrac- 
tion envelope and the scale factor appearing in Eqs. (4) 
and (5). This means that for an arbitrary state p, if it is 
possible to postselect it in the same state by measuring 
in different points {x,z) and (x',z'), in general we will 
have for the detection probability density [Eq. (7)] 



s(x,z) = \ip{x,z)){ip{x,z)\, 



(6) 



g{x,z) ^ g{x',z'). 



(8) 
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FIG. 3: (Color online). Bloch sphere (in the Hammer- Aitoff 
map projection) for the spatial qubits. Each curve corre- 
sponds to the trajectory traced out by the state vectors one 
can measure by displacing the detector in the x direction for 
the detection planes z shown in the legend. The latitude 
6 € [0, tt] and the longitude (p G [— vr, n] parametrizes an ar- 
bitrary pure state IV-} = cos(6'/2)|0 + e'*^ sin(6'/2)|r) on the 
sphere surface. 



2. Expected results when Alice uses a point detector 

Let us go back to the RSP protocol of spatial qubits 
when Alice uses a point detector in the setup of Fig. 2(a). 
If photon A, which is part of the entangled pair in the 
state of Eq. (1), is detected at the position {x,z), then, 
from Eqs. (3) and (6), the state of Bob's spatial qubit 
will be given by 



oc ifi*{x,z)\l) + ipi{x,z)\r) 



(9) 



after she has transmitted 1 cbit to him, namely the out- 
come of her measurement. Accordingly, if Alice wants to 
prepare Bob's qubit in the pure state of Eq. (2), then she 
must place her detector in the position (x, z) such that 



ipi{x, z) (X l3* and ipr{x,z)(xa* 



(10) 



As discussed before, with this procedure she can measure 
arbitrary pure states so that she can remotely prepare ar- 
bitrary pure states as well. After receiving the cbit from 
Alice, Bob does not need to perform a unitary correc- 
tion since, with this strategy, the target state has been 
already prepared as it was intended to be. From Eq. (7), 
the probability density for this preparation will be 



g{x,z) = TTAB{[s{x,z)^iB]\'f)AB{'^\} 
\ipr{x,z)\'^ + \ipi{x,z)\'^ 



(11) 



where Is is the identity operator for Bob's qubit, and 
ifj{x,z) (for j — l,r) is given either by Eq. (4) or (5) 
depending whether Alice measures in region II or I of 
Fig. 2, respectively. It is worth mentioning that in this 
idealized situation, the fidelity and purity of the remotely 
prepared spatial qubits would be both equal to one. 



B. Realistic measurement: Detector with finite 
resolution 

In the real world, any detector has a finite spatial res- 
olution, which means that the best it can do is to locate 
the photon within some narrow interval Ax around the 
position X. This single-pixel detector is illustrated in the 
inset of Fig. 2(a). For this device, the associated mea- 
surement operators can no longer be written as in Eq. (6). 
A realistic measurement operator can be derived by not- 
ing that the probability of detecting the photon within 
the range {x — Ax/2,x + Ax/2) in a given transverse 
plane z is 



Pix,Ax,z) 



x+Ax/2 



g{x' , z) dx' . 



(12) 



t:-Ax/2 

Hence, using Eq. (7) is straightforward to show that 



P{x,Ax,z) — Tt S{x,Ax,z)p 



where 



S(x, Ax, z) 



x+Ax/2 



-Ax/2 



s{x' , z) dx' , 



(13) 



(14) 



is the measurement operator for a detector with finite 
resolution [36]. 

Let us see what are the consequences for the RSP pro- 
tocol of spatial qubits, when Alice uses this realistic de- 
tector in her apparatus. In the same way as before, if she 
wants to prepare Bob's qubit in the target state given by 
Eq. (2), she must place her detector at the position {x, z) 
satisfying the condition (10). The two-photon state in 
Eq. (1) after Alice's measurement is transformed as fol- 
lows: 



1^' 



AB 



S{x, Ax, z) 



1/2 



1; 



AB 



AB{'^\[Six,Ax,z) 



(15) 



IbII*) 



AB 



where the existence of 



S{x, Ax, z) 



1/2 



is guaranteed by 



the positivity of S{x, Ax, z) [36]. Once she communicates 
to Bob the outcome of her measurement, his spatial qubit 
is left in the state 



PB = TrA{W)AB{-^'\) 

^rr{x. Ax, z) 
^lr{x,Ax,z) 



N 



^rl{x, Ax, z) 



(16) 



where N — [$,.^(2;, Ax, z) + ^ii{x. Ax, z)]~^ , is a normal- 
ization constant and the matrix elements Ax, z) 
are given by [using Eqs. (3), (6), and (14) ] 

^ij{x,Ax,z) = {i\S{x,Ax,z)\j) 



x-Ax/2 



x-Ax/2 



(Pi{x', z)ip* {x, z) dx', (17) 
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for i,j = r, I. This representation will be useful to avoid, 
as much as possible, cumbersome expressions. 

As intuitively expected, a detector with finite resolu- 
tion introduces impurity in the remotely prepared spa- 
tial qubit because of the imperfect filtering it provides. 
Mathematically, this can be seen from Eq. (16) where the 
functions ^ij{x, Ax, z) defined in (17) cannot factorize in 
a product of functions like (pr{x, z) and ^pi{x, z). This will 
only happen in the limiting case of a point detector [25] 
where we have 



iPi{x,z)ip*{x,z). 



(18) 



Despite this, the above condition can be approximated as 
long as the pixel detector width be small compared with 
the scale where the functions like (4) or (5) oscillate [25] . 
At the other extreme, there is a detector with infinite 
width. In this case, thanks to the orthogonality of the 
probability amplitudes ipi{x,z) [34] in Eq. (4), one can 
show that 



lim $ij(x. Ax, z) 

Aa;— J-CJO 



(19) 



In this case, Alice would remotely prepare only the max- 
imally mixed state pb = 5 Is- This is also an expected 
result, since a completely opened detector, or a "bucket" 
detector, does not register where the photon has arrived 
and therefore is equivalent to trace over Alice's photon 
in Eq. (1) [22]. 



1. Figures of merit of the protocol 

To evaluate the performance of the RSP of spatial 
qubits according to Alice's measurement strategy, we 
have chosen three figures of merit, namely, the proba- 
bility of preparation, the fidelity and the purity of the 
remotely prepared states. 

The probability of preparation is obtained by replacing 
Eq. (11) into (12). Using the definition of Eq. (17) it can 
be written as 

P(x,Ax,z)^ ''-^^'^^'"^ + ''"(^'^^'"^ (20) 

To calculate the fidelity between the state remotely 
prepared by a realistic detector [pB in Eq. (16)] and the 
target state remotely prepared by an idealized point de- 
tector [[?/;) B in Eq- (9)], we must normalize the latter and 
then, using the definition F = ('0|/5|V'), we get 



F(x, Ax, z) — 



E»j^r,i y^iei (a^: {x, Ax, z)ip*^^{x, z) 

J2mn=r I I'Pmix, z)|2$„„(x, Ax, z) 

(21) 

where © indicates addition modulo two (l (B 1 — r and 
r© 1 = 0- 

Finally, the purity of the remotely prepared qubits in 
Eq. (16), using the definition V = Tt{p^), is given by 

[X, AX, z) - j^^^^^^ ^ ^^^^^^ ^^^^ . [ZZ) 



With the fidelity and purity, one can quantify the ef- 
fects of the finite resolution of the detection system over 
the remotely prepared states. It is easy to check from 
Eqs. (21) and (22), respectively, that both quantities are 
smaller than one for a detector of finite width Ax. The 
only exception is for the remote preparation of qubits in 
the logic basis {[/), |r)} when Alice measures at the im- 
age plane with a detector whose width. Ax, is not large 
enough to "see" both slits simultaneously. In this case, 
the fidelity of the remotely prepared state and its purity 
will be one. This happens because when one integrates 
over the detector width there will be only a change in the 
phase which is global and thereby, irrelevant. 

In the limiting case of a point detector when Eq. (18) 
holds, the fidelity in Eq. (21) and the purity in Eq. (22) 
goes to one for every remotely prepared state, whereas 
in the case of a bucket detector, the condition (19) is 
satisfied and these quantities will be both 1/2. 



2. Comparing the expected fidelities with the experimental 
results by Taguchi et al. 

We can compare the theoretical values for the fidelity 
of the remotely prepared spatial qubits with those ob- 
tained in the experiment performed by Taguchi et al. in 
Ref. [17]. In that work, the authors have demonstrated 
the remote preparation of six states of spatial qubits, 
through spatial postselection as described in this section. 
The parameters characterizing their target states [a and 
/3 of Eq. (2)] are shown in the first and second columns 
of Table II. They were obtained from Alice's detector po- 
sitions (x, z) provided in [17] and the parameters a, d, A 
and / they used (shown in the caption of Table II) . Re- 
placing these values in Eqs. (3) and (4) and normalizing 
the states, we can determine a and /3. Similarly, taking 
into account the detector width in [17] (Ax — 20 /im) 
and using Eqs. (4), (17), and (21), we calculated the ex- 
pected values of the fidelity between these target states 
and the remotely prepared ones. The obtained values are 
shown in the third column of Table II, while the fourth 
column shows the experimental results of [17]. One ob- 
serves clearly that the experimental values are within 
the bound imposed by the theoretical ones and further- 
more, following their tendency: the states with expected 
higher fidelities showed higher fidelities (first and sixth 
rows), those with expected intermediate values showed 
intermediate fidelities (second and fifth rows), and the 
same for the states with expected lower fidelities (third 
and fourth rows). The reduction of the experimentally 
obtained fidelities, in comparison with the theoretically 
expected, can be attributed to several factors: the non- 
perfect preparation of the maximally entangled photon 
pair, the error in the positioning of Alice's detector and 
the tomographic process of Bob which also includes a de- 
tector with finite resolution. However, even these sources 
of errors were negligible in [17], the bound on the maxi- 
mum achievable fidelities would be set by the finite reso- 
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TABLE II: Comparison between the theoretical fidelities cal- 
culated from Eq. (21) and the experimental ones obtained in 
Ref. [17]. The first and second columns show the parameters 
which characterizes the target states [see Eq. (2)] intended 
to be prepared in [17]. The experimental parameters of [17] 
are: a = 20 ^m, d — 150 fim, X — 810 nm, f — 5 cm, and 
Aa; = 20 fim. 

Target state |?/))b Fidelity 

|q| arg(^/Q) Theory Exp. [17] 



0.979 


-0.113 


98.4% 


88.7% 


0.776 


-2.319 


89.9% 


86.1% 


0.742 


-1.159 


89.6% 


84.1% 


0.670 


1.159 


89.6% 


84.1% 


0.631 


2.319 


89.9% 


87.1% 


0.201 


0.113 


98.4% 


91.2% 



lution of Alice's detector within the optical system. 

C. Numerical analysis 

We have analyzed, numerically, the figures of merit of 
the RSP protocol described so far (see Fig. 2) and here 
we present and discuss our results. The parameters used 
in our simulations are shown in Table I and their val- 
ues are typically employed in experiments. In particular, 
the value we have chosen for the detector width, together 
with the other parameters, allows for a very good spatial 
filtering without compromising too much the detection 
rate [27, 28]. Our simulations were performed in the fol- 
lowing steps: 

1. The surface of the Bloch sphere is divided into a 
grid of 300 x 600 pixels {9i x which give us 
1.8 X 10^ states to be analyzed. 

2. We choose one figure of merit to maximize. Here 
this will be either the probability [Eq. (20)] or the 
fidelity [Eq. (21)]. 

3. Alice's detector is scanned in the x direction 
in steps of 20 /im (= Ax) within the range 
[— 7r///ca, 7r//fca] corresponding to the first-order 
diffraction zone at the focal plane of the lens. Thus 
for each transverse plane we get 250 transverse po- 
sitions Xm (m = 1, . . . , 250). 

4. This scanning is performed along 3 x 10^ transverse 
planes, z„ {n = 1, . . . , 3 x 10^), within the range 
[/, 2/] , which gives us a resolution of 1 /im in the 
longitudinal direction. 

5. For a given plane z„, each point (x^, z„) defines an 
(idealized) remotely prepared state given by Eq. (9) 
and after normalization we determine its corre- 
sponding pixel {9i^(j)j) in the Bloch sphere. The 
figure of merit of this state is then calculated and 
the obtained value is inserted into the pixel {6i, (pj). 



TABLE III: Comparison between the average, minimum, 
maximum values and standard deviation of the probability 
(P), fidelity (_F) and purity [V). The second (third) column 
is associated with RSP via spatial postselection when Alice 
maximizes the probability (fidelity) and was obtained from 
Figs. 4(a)-4(c) [Figs. 4(d)-4(f)]. The fourth column shows 
the statistics for RSP via generalized measurement (Sec. IV) 
obtained from the results shown in Figs. 6(a)-6(c). 





Spatial postselection 
Max. probability Max. fidelity 


POVM 


{P) 


0.02952 


0.01871 


0.01547 


p . 

mm 


0.0001 


0.0001 


0.01521 


P 

-1 max 


0.1978 


0.1978 


0.01559 


ap 


0.04728 


0.03310 


0.0001 


(F) 


0.8764 


0.9024 


0.9997 


-fmin 


0.1964 


0.2046 


0.9995 


-^max 


1 


1 


1 


ap 


0.1568 


0.1335 


0.0002 


(V) 


0.8853 


0.8987 


0.9995 


Pmin 


0.5057 


0.5058 


0.9989 


Pmax 


1 


1 


1 


a-p 


0.1273 


0.1183 


0.0004 



6. In the next plane Zn+i the procedure is repeated. 
The blank pixels are occupied with the calculated 
values of the figure of merit. If some pixel was al- 
ready occupied from a scan on previous plane (s), 
then the figure of merit with higher value will over- 
ride the smaller one. Thus, for a given vector which 
can be measured in different points, only the one 
with the highest value of the chosen figure of merit 
will remain. 

This procedure is repeated for all planes z„ so that in the 
end we are left with the highest values for the figure of 
merit we intended to maximize. The two remaining fig- 
ures of merit are calculated at the points where the first 
one has been maximized. It is important to mention that 
for f Zn < 1.8/, that is, in the Fraunhofer approxima- 
tion [see Fig. 2(b)], both the state (9) and the figures of 
merit are calculated with the help of the probability am- 
plitude given by Eq. (5), otherwise, for 1.8/ < Zn < 2/ 
we used Eq. (4). 

Firstly, we searched for the higher probabilities of re- 
mote preparation and the obtained results are shown in 
Fig. 4(a). The corresponding fidelities and purities of 
the remotely prepared states are shown in Figs. 4(b) and 
4(c), respectively. Secondly, we searched for the higher fi- 
delities of the remotely prepared states and the obtained 
results are shown in Fig. 4(d). Figures 4(e) and 4(f) 
show the corresponding probabilities and purities, respec- 
tively. Table III shows the statistics for the three figures 
of merit taken from the graphics in Figs. 4(a)-4(c) (sec- 
ond column) and Figs. 4(d)-4(f) (third column). 
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FIG. 4: (Color online). Numerical results for the three figures of merit of RSP of spatial qubits through spatial postselection. 
The first row refers to Alice's choice of maximizing the probability of preparation and the result is shown in the Bloch sphere 
(a), (b) and (c) show the corresponding fidelity and purity of the remotely prepared states, respectively. In the second row, 
Alice maximizes the fidelity of the prepared states (d) and their corresponding probability of preparation and purity are shown 
in (e) and (f), respectively. The procedure for the simulation is described in Sec. IIIC. 



1. Discussion 

The effects of diffraction and the scale factor of its 
associated pattern are clearly observed by the nonuni- 
form probability distribution in either case [Figs. 4(a) 
and 4(e)]. Also, the effect of a detector with finite spa- 
tial resolution can be seen by the non-unity values for the 
fidelity [Figs. 4(b) and 4(d)] and purity [Figs. 4(c) and 
4(f)] of the remotely prepared states. 

When Alice looks for the maximal probabilities, the 
higher values are found in the poles of the Bloch sphere 
and at their vicinity. This means that the correspond- 
ing states are measured at the image plane and planes 
close to that, where the hght intensity is stronger. For 
the majority of the superpositions between \l) and |r), 
however, the higher probabilities are found far from the 
image plane so that the achieved values are strongly re- 
duced by diffraction, as can be seen by the blue zones in 
Fig. 4(a). Quantitatively, these results are show in the 
second column of Table III. 

For several states in the Bloch sphere, the detector po- 
sition (x, z) where the maximal probabilities were found, 
are not the same where the fidelity is maximal. This 
can be seen, by comparing the graphics in the first row 
with those of the second row of Fig. 4 or, comparing the 
second and third column of Table III. In the positions 
where Alice finds the higher fidelities, the corresponding 
probabilities of preparation decrease while the purities 
increase. These results are so because in order to maxi- 
mize the fidelities, she must place her detector in points 
where the spatial filtering is stronger, which affects both 



the probability and purity as described above. The max- 
imum value of one achieved by the fidelity and purity, 
whatever the figure of merit Alice wants to maximize, is 
obtained only when she remotely prepares \l) or \r) as we 
discussed in Sec. IIIB. 

In our analysis we have considered a detector with fixed 
width. It would be possible for Alice to use a detector 
with variable width. This would allow her to improve 
the fidelity (or purity) of the remotely prepared states. 
For instance, let us consider the case where she has max- 
imized the fidelities shown in Fig. 4(d). It can be seen 
there that the lower fidelities are mainly distributed over 
the regions between the poles and the equator of the 
Bloch sphere. At the same time, in those regions the 
probabilities of preparation are, in general, very small. 
Therefore, reducing the detector width to increase the 
lower fidelities would, simultaneously, reduce the already 
low probabilities to prohibitively low levels. 



IV. RSP OF SPATIAL QUBITS VIA 
GENERALIZED MEASUREMENTS 

In this section we will present an alternative measure- 
ment strategy with which Alice can remotely prepare 
Bob's spatial qubit. This strategy is based on our re- 
cent work [30] and, as we will see, at the expense of in- 
creasing the classical communication cost as well as the 
experimental resources, it improves the figures of merit 
of the RSP protocol when compared with spatial postse- 
lection alone. The strategy can be divided into two steps: 
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FIG. 5: (Color online). Schematic representation of the mea- 
surement strategy for RSP of spatial qubits proposed in this 
work. In the dashed box, a half-wave plate (HWP) behind 
each slit followed by a polarizing beam-splitter (PBS) imple- 
ments the two-outcome POVM given by Eq. (23) [30]. In each 
output p, a two-pixel detector is placed at the focal plane of a 
lens and it is allowed to displace in the x direction. For appro- 
priate choice of the HWP angles [Eq. (35)] and the detector 
positions [Eqs. (30), (39) and (40)], the remotely prepared 
state is given either by (41)-(44) for point detectors or by 
(46) for detectors with finite resolution. The left side inset 
shows a sketch of the two-pixel detector and the fixed separa- 
tion between the pixels. The right side inset shows the range 
where the two-pixel detector moves (the shaded region). The 
dashed line shows how the phase changes according with the 
detector position. 



first, the implementation of a two-outcome POVM with 
the polarization as the ancillary system and second, the 
spatial postselection at the focal plane of the lens by a 
two-pixel detector in each output of the POVM. In the 
following we describe in detail the whole process. 



A. The POVM and its physical implementation 

In the first step of the new strategy, Alice will imple- 
ment a two-outcome POVM on her photon. The physical 
implementation of a POVM requires the extension of the 
Hilbert space of the system to be measured, which can 
be provided by an ancillary quantum system, or ancilla 
[8]. Following our proposal [30], the ancilla for Alice's 
spatial qubit will be the photon polarization. For sim- 
plicity, let us assume that the overall two-photon state 
is |*)ab «) \H)a\H)b, where is given by Eq. (1) 

and \H)i {i = A, B) denotes horizontal polarization. The 
dashed box in Fig. 5 sketches the setup for implement- 
ing the POVM. Conditional to the passage of photon A 
through the slit j, its polarization is rotated by 29 j with 
a half- wave plate (HWP). This unitary operation entan- 
gles the polarization and spatial degrees of freedom. Af- 
ter that, a polarizing beam-splitter (PBS) measures the 
polarization in the basis {\H) , \ V)} (where V represents 
vertical polarization), accomplishing the POVM on the 



spatial qubit. 

Let lip be the elements of the POVM described above, 
where the subscript p — H,V labels the two possible 
outcomes. Mathematically, they can be written as 



A^pAp 



{H\U^\p){p\U\H). 



(23) 



One can easily check that all properties of a POVM (her- 
miticity, positivity and completeness) are satisfied. The 
unitary operator is given by [37] 



j=r,l 



where 



m) = 



cos(26'j) 
sin(26lj) 



-sin(26'j) 
cos(26lj) 



(24) 



(25) 



describes the conditional polarization rotation. There- 
fore, when Alice and Bob share the two-spatial qubit 
state in Eq. (1) and she implements the POVM (23) on 
her photon, the two-photon state if the outcome is p will 
be transformed as follows: 



1 



>AB 



{Ap' 



1} 



I*) 



AB 



2^/K 



{{p\ei)\l)A\r)B 



{pmr)A\l)B), 

(26) 



where 



|0„) =cos(20„)|i/)+sin(20„)|l^), 



for n — l,r and 



Pr, 



\{pm' + \(p\Or)\' 



(27) 



(28) 



is the probability that photon A exits through the output 
p, with Ph + Pv = 1. 

The unitary operation of Ec^. (24) may be experimen- 
tally implemented with the help of programmable liquid- 
crystal displays (LCDs) instead of HWPs behind the slits 
(Fig. 5). LCDs are multipixel optical devices which can, 
in a controlled way, rotate the light polarization accord- 
ing to the pixel's variable phase retardance. Recently, 
they have been used to manipulate the amplitude [29, 38] 
and/or phase [27, 28] of spatial qudits with high accuracy. 



B. Measurement at the focal plane by two-pixel 
point detectors 

1. The two-pixel detector 

After being transmitted or reflected at the PBS, pho- 
ton A propagates toward the detection system which is 
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placed at the focal plane of a lens, as shown in Fig. 5. 
Likewise Sec. Ill, here we will consider first the ideal- 
ized situation where Alice has point detectors in her 
setup, and after that we discuss the realistic case of de- 
tectors with finite resolution. Now, rather than a single 
pixel detector, Alice's detection system (in each output 
p = H,V) will be comprised of two pixels (labeled by 
1 and 2) which are mounted over the same platform, so 
that they can move simultaneously in the transverse di- 
rection a;. In addition, we will impose that the pixels 
are separated by a fixed transverse distance in order to 
postselect orthogonal states. Using Eqs. (3) and (5) it is 
easy to show that this condition is satisfied when 

{(p{xip,f)\Lp{x2p,f)} oc sine ^ ^"J^^ ^ sine ^ ^"'^'^p 

'k {x2p - xip)d 
2/ 



0, 



(29) 



where Xjp for j — 1,2 denotes the transverse position 
of detector j at the output p. On one hand it will be 
satisfied at the zeros of the sine functions, which are as- 
sociated with the zeros of the diffraction envelope, and 
thus does not concern us. On the other hand, from the 
cosine function (associated with the interference pattern) 
we can obtain the fixed distance between the detectors 
satisfying (29), which is X2p — xip = irf /kd as shown 
in the left inset of Fig. 5. Due to the periodicity of the 
cosine function we can further constrain the two-pixel 
detector to move in a single period within the interval 
[—TTf/kd,Trf/kd] as shown in shaded region of the right 
inset of Fig. 5. As can be seen in this inset, the detection 
efficiency due to the diffraction envelope is higher in that 
region. From this discussion, the relation between the 
detector positions in each arm p can be written as 



X2p 



':ip ' 



kd' 



(30) 



where ® denotes that the addition is performed mod- 
ulo 27r/ / kd, ensuring that both detectors are confined to 
move in the interval [—TTf/kd,TTf/kd], keeping the sepa- 
ration of nf/kd as shown in insets of Fig. 5. 

The measurement operator associated with the point 
detector j in the focal plane of the lens and at the output 
p of the PBS (see Fig. 5) is given by s{xjp, /), as defined 
by Eq. (6). Using Eqs. (3) and (5), the corresponding 
postselected (non-normalized) state is 

Mx.pJ)} = e^^'^'^^/'^M^jpJ) (lO +e'''""/^|r) 

(31) 

Therefore, the two-photon state in Eq. (26) after Alice's 
measurement will be transformed as follows: 



\'^jp)aB oc \ip{XjpJ))A\ll^jp)B, 



(32) 



is the remotely prepared state of Bob's photon. The vec- 
tors \9n) (for n = r,l) are given by Eq. (27). Taking into 
account the probability for the outcome p [Eq. (28)] in 
the POVM, the probability density for this preparation 
will be 

f) = PpXTrABm,p)AB{^jp\} 

\{pm' + \{p\9r)\' , .,,2 

= ^ m{xop,j)\ ■ (34) 

Equation (33) makes clear the roles of the POVM and 
the detector position for the realization of RSP of spatial 
qubits. The former will set the real part of the ampli- 
tudes a and /3 of the target state (2) while the latter, the 
relative phase arg(/3/Q!). This will be discussed in more 
detail next. 



2. Conditions for remotely preparing the target state 

Having described the POVM and the detection system, 
we now turn to the conditions that Alice has to accom- 
plish to remotely prepare Bob's photon in the target state 
given by Eq. (2). Firstly, she must set the wave plates 
behind the double slit as follows: 



9/ = e = - cos" 
2 



^(l«l), 



(35) 



Thus, using Eqs. (27) and (33), the photon B will be 
prepared in the state 



\'^]h)b 

\^3v)b 



sin(2e)|0 +e''=''"^«//cos(2e)|r), (36) 
cos(2e)|0 + e'^'^^^^'i sin(2e)|r), (37) 



given that photon A is measured by the detector j = 1,2 
at the exit port p ~ H,V oi the PBS. In addition, the 
probability density for the preparation (34) becomes 



eixjpj) = 



I'piixjpJW 



(38) 



The relative phase of the target state, e.rg{f3/a), is 
controlled by the detector positions. Previously, we saw 
how the two-pixel detectors move within each arm [see 
Eq. (30)]. Now, we impose an extra condition between 
the positions of these arms {p = H, V) given by 



XjH 



-XjV, 



(39) 



such that the whole detection system becomes con- 
strained by Eqs. (30) and (39). Therefore, the choice 
of the position of one detector will immediately set the 
position of the remaining three. Thus, if 



where 



XlV 



kd 



arg 



(40) 



\iPjp)b oc {p\er)\l) + e^^'^'-^'"/^ {p\9i)\r), (33) we get xih = -xiv and X2v = xiv ® T:f/kd ■ 



-X2H- 
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Accomplishing the conditions (35) and (40) , Alice will, 
conditioned to detect her photon with detector j = 1,2 
at output p — H,V, remotely prepare Bob's spatial qubit 
as follows: 

a\l)+P\r)= (41) 
\Mb^ a\l)-P\r)= <7,|V)b, (42) 
\^ih)b^ m+a\r)= ct^Wb, (43) 
\Mb^ m-a\r)^ cTy\ib)B, (44) 

where (for n = x,y,z) denotes the Pauli operator 
which Bob must apply on his qubit after receiving the two 
bits of classical communication from Alice, in order to 
recover the target state [Eq. (2)] she intended to prepare. 
It is easy to see that with this strategy Alice can prepare 
arbitrary pure states. First, the POVM which sets the 
real parts of the amplitudes, allows for values |a| € [0, 1]. 
Second, as shown by the dashed line in the right inset 
of Fig. 5, the range where the detectors move suffices for 
imparting relative phases from — tt to tt. As in the case 
of spatial postselection studied in Sec. Ill A, the fidelity 
and purity of the remotely prepared qubits would be both 
equal to one, if the detectors were just points. 



Measurement at the focal plane by two-pixel 
detectors with finite resolution 



Let us now turn to the realistic case where Alice's de- 
tection system is composed by detectors with finite res- 
olution. The whole procedure is exactly as we just de- 
scribed and can be summarized as follows: 



1. Alice implements the two-outcome POVM (23) on 
her half of the maximally entangled photon pair, 
setting the wave plates satisfying conditions (35). 

2. Alice fixes her detection system, in each output p, 
at the focal plane the lens and satisfying conditions 
(30), (39), and (40). 

3. Alice communicates 2 bits of classical information 
to Bob, conditioned to detect her photon with de- 
tector j = 1, 2 at output p = H,V. 

4. Bob performs the appropriate unitary transforma- 
tion given by (41)-(44) according with the result 
communicated by Alice. 

The two-photon state (1) after Alice has implemented 
the POVM (23) is transformed in (26), whether the out- 
come was p. Thus when she measures photon A with the 
detector with finite width Ax, this state is transformed 
as follows: 



S{xjp,Ax,f) 



1/2 



«) li 



AB, 



(45) 



where S{xjp, Ax, /) is the measurement operator, defined 
in Eq. (14), for a reahstic detector j at the output p. Once 
Alice communicates to Bob the outcome of her measure- 
ment, his spatial qubit is left in the state determined by 
fP^ = 1:,a{W„)ab{^'J)- Using Eqs. (14), (17), (27) 
and (35) one can show that 



Pb ■^ 



^ sin2(2e) <I>„(a;jH, Ax,/) sin(2e) cos(2e) $w(a;jff , Ax, /)^ 

^sin(2e)cos(2e)$;,(x,j^,Ax,/) cos2(2e) $„(x,ff , Ax, /) ^ 

^ cos2(2e)'I>„(xjy, Ax,/) s\n{2Q)cos{2Q)^ri{xjv,AxJ)\ 

sin(2e)cos(2e)$i^(xjy,Ax,/) sin2(2e) $h(xjv. Ax, /) ) 



if p^H, 



if p^V. 



(46) 



For point detectors, when the Eq. (18) holds, p-'g reduces 
to \tpjp)B in Eqs. (41)-(44). For bucket detectors the 
condition (19) holds and Alice would be able to prepare 
only pure states lying in the poles of the Bloch sphere 
{|/), |r)}, or mixed states along the straight line connect- 
ing these poles. This already improves over the spatial 
postselection in which a bucket detector would prepare 
just maximally mixed states (see Sec. IIIB). 



D. Numerical analysis 



Similarly to Sec. Ill, here we have also analyzed, nu- 
merically, the figures of merit of the RSP protocol when it 
is performed via generalized measurements (see Fig. 5). 
The parameters used in our simulations, shown in Ta- 
ble I, are the same as those used in the previous section. 

The first figure is the probability of preparation. This 
is obtained first by integrating the probability density 
(38) over the detector width Ax. Using the definition 
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FIG. 6: (Color online). Figures of merit for RSP of spatial 
qubits via generalized measurements (see Fig. 5). (a) Total 
probability of preparation; (b) fidelity and (c) purity of the 
remotely prepared states. 



(17) we get 



P{xjp,Ax,f) 



(47) 



Since now Alice has four detectors which, upon a click, 
always prepare Bob's qubit up to a unitary correction, 
the total probability of preparation will be 



^tot(a;iy, 

j = l p=H,V 



2a 



X1V + 



Ax 



ka 



dx' 



sine 



/ kax' 



sine I y (2^' ® T^f/kd) 



(48) 



Note that the total probability is written as a function 
of only one detector position, xiy, since this sets the 
position of the others through (30), (39), and (40). Fig- 
ure 6(a) shows the plot of this probability which varies 
around 1% due to the diffraction envelope (see the right 
inset of Fig. 5). 

The fidelity of the remotely prepared states must be 
calculated between the target state prepared when a 
point detector j = 1, 2 clicks at output p = H,V [|V'jp)-B 
in Eq. (33)] and the state effectively prepared with real- 
istic detectors [fPg in Eq. (46)]. Doing so we get 



F{xjp, Ax, 9) = cos4(2e) -t- sin''(2e) 



2cos2(2e)sin^(2e) 
^n{xjp,Ax,f) 



Xj^p-\- Ax/2 



-Ax/2 



dx' 1-^/(2;', /)p cos 



/fcd(a 



(49) 



The purity of the remotely prepared states in Eq. (46) 
will be given by 



Vixjp,Ax,e) = cos^(2e) 
1 



sin^(4e) 



sin''(2e) 

^IriXjp, 



^u{xjp,Ax, f) 



(50) 



It turns out that the variation of both quantities, fidelity 
and purity, within the range where the de- 

tectors are allowed to move is negligible. This can be 
verified by replacing the parameters of Table I into the 
terms of Eqs. (49) and (50) which depend on Xjp and 
calculating those terms in that range. Therefore the fi- 
delity and purity of the remotely prepared states can be 
written, respectively, as 

F(e) = 1 - 0.002[cos(2e) sin(2e)]^ (51) 



V{Q) = 1 - O.G04[cos(2e) sin(2e)]^ 



(52) 



and their plots as a function of the HWP angles given by 
Eq. (35) are shown in Figs. 6(b) and 6(c). The behav- 
ior of both curves can be understood as follows. When 
e = ±7r/4 or e = each output of the POVM com- 
pletely specify from which slit the photon came. In this 
case, the states are prepared with unity fidelity and pu- 
rity no matter the detector width, similarly as the mea- 
surement at the image plane in spatial postselection (see 
Sec. III). On the other hand, for 8 = ±7r/2 there is 
no "which-slit" information whatsoever at the outputs of 
the POVM. Thus, the fidehty and purity depends only 
on the detector width, reaching their minimum values. 
The constant term of Eqs. (51) and (52) accounts for 
this dependence. For distinct values of O there is par- 
tial which-slit information and interference so that the 
fidelity and purity took intermediate values. 
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The statistics for the three figures of merit obtained 
from the results of Figs. 6(a)-(c) is shown in the fourth 
column of Table III. As it is seen in the graphics of Fig. 6, 
the variation of the three quantities is very small and, 
accordingly, their standard deviation is negligible. In the 
next section we discuss these results by comparing them 
with those obtained from spatial postselection. 

V. COMPARISON OF THE METHODS 

In this section we compare the two measurement 
strategies for RSP of spatial qubits discussed here, re- 
garding the required resources for each one as well as 
their figures of merit. 

A. Classical communication cost and experimental 
resources 

Our proposed strategy discussed in Sec. IV requires 1 
ebit and 2 cbits to accomplish the RSP and thereby, in- 
crease the classical communication cost of the protocol 
comparing with spatial postselection (Sec. Ill) which re- 
quires 1 ebit and 1 cbit. Resorting to the comparison 
between RSP of spatial and polarization qubits estab- 
lished in the Introduction, this is expected since spatial 
postselection is the analog of a probabilistic implementa- 
tion in polarization while our strategy is the analog of a 
deterministic one. 

In terms of experimental resources, our generalized 
measurement strategy is much more expensive than spa- 
tial postselection, in the same way that deterministic 
RSP of polarization qubits [16] is more expensive than 
probabilistic implementations [13]. For spatial postselec- 
tion, Alice needs a single-pixel detector moving between 
the focal and image plane, and a logic circuit to send the 
classical message to Bob. On the other hand, with our 
proposed strategy Alice needs an apparatus for perform- 
ing the POVM (which can be done by means of a LCD 
[28]), four detectors and a more complex logic circuitry 
to encode the 2 cbits which will trigger Bob's unitary 
correction [see Eqs. (41)-(44)]. This last step may be the 
more challenging in a possible experimental implementa- 
tion but it is doable in principle. 

B. Figures of merit 

1. Probability 

In Table III it is seen that the average values for the 
probability of preparation are higher when Alice employs 
spatial postselection. However, looking for the maximum 
and minimum values of this probability and the standard 
deviation of its distribution over the Bloch sphere one can 
ask how many states are in fact prepared with higher 
probabilities comparing with our proposed generalized 




FIG. 7: (Color online). Bloch sphere comparing the prob- 
abilities for RSP according to Alice's measurement strategy. 
The red color indicates the states which are prepared with 
higher probabilities when Alice adopts spatial postselection. 
In blank, the states prepared with higher probabilities when 
Alice adopts the generalized measurement. The comparison 
in (a) refers to Alice's choice of maximizing the probability 
of preparation [Fig. 4(a)] and (b) the fidelity of the remotely 
prepared states [Fig. 4(e)]. 



measurement. This comparison is shown in Fig. 7, where 
the red color indicates the states prepared with higher 
probabilities by spatial postselection, and the blank in- 
dicates the states prepared with higher probabilities by 
generalized measurements. In Fig. 7(a) the comparison 
refers to Alice's choice of maximizing the probability of 
preparation via spatial postselection [see Fig. 4(a)]. In 
this case, for 61.4% of the states analyzed, the general- 
ized measurement strategy outperforms spatial postselec- 
tion. In Fig. 7(b) the comparison refers to Alice's choice 
of maximizing the fidelity of the remotely prepared states 
via spatial postselection [the corresponding probability is 
shown in Fig. 4(e)]. In this case, 70.2% of the states an- 
alyzed presented higher probabilities when prepared via 
generalized measurements. Therefore, our measurement 
strategy provides higher probabilities of remote prepara- 
tion for the majority of the states when compared with 
spatial postselection, whatever the figure of merit Alice 
wants to maximize. 
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2. Fidelity and purity 

The fidelity and purity of the remotely prepared states 
are the same for both measurement strategies only when 
Alice prepares states lying either in the poles or in the 
equator of the Bloch sphere. Otherwise, the generalized 
measurement strategy always provides higher values for 
these figures of merit. In fact, in our proposed strategy 
all the states are prepared with fidelity and purity close 
to the ideal value of one (achieving one for the prepara- 
tion of |/) or |r)), as can be seen in Figs. 6(b) and 6(c), 
respectively. On the other hand, for spatial postselec- 
tion many states are prepared with fidelity and purity 
much lower than the minimum values achieved with gen- 
eralized measurement. Quantitatively, this can be seen 
by comparing the statistics for these quantities shown in 
Table III. 

As discussed in Sec. Ill C, when adopting spatial post- 
selection, Alice could employ a detector with variable 
width to increase the fidelities (or purities) in the posi- 
tions which they are low. Doing so, she could, in princi- 
ple, obtain results as good as those obtained by general- 
ized measurement. However, the required spatial filter- 
ing would be so strong that the probability of preparation 
would become inviable. 



VI. RSP OF SPATIAL QUDITS 

Let us briefly discuss the RSP protocol for spatial qu- 
dits. In Ref. [25], Taguchi et al. have experimentally 
demonstrated the remote preparation of spatial qutrit 
states using spatial postselection. In this section wc dis- 
cuss the limitations of this method for RSP arising from 
the increasing of the qudit dimension and, accordingly, 
the number of free parameters of the state to be prepared. 
Additionally, we present a generalized measurement [30] 
which will overcome these limitations. 

To start with, it is well known that a SPDC source 
can generate maximally entangled states of spatial qudits 
[20, 21]. The setup for this is similar to that shown in 
Fig. 1 with the double slits replaced by D-slits arrays. 
When the pump beam is focused at the plane of these 
apertures {zap) the two-photon state after transmission 
through them will be [20, 21] 

l*)^^ = 77! E ^'''\j)a\-j)b, (53) 

where i = (I? — l)/2 and — kdP p /2zi^p is a phase fac- 
tor that depends on fixed experimental parameters and 
in principle may be compensated by introducing external 
phase shifts into the slits. Let us assume that this is the 
case so that in the following we consider jjLj — OVj. 

Now, Alice's goal is to perform a measurement on her 
half of the maximally entangled pair, communicate the 
outcome to Bob, and then remotely prepare his spatial 




^ HWP [ Phase shifter 



FIG. 8: (Color online). Schematic representation of the mea- 
surement strategy for RSP of spatial qudits proposed here. 

qudit in the pure state 

e 

- E ^^l^')' (54) 

with X]^=-f|cjP = 1- Let us again consider the two 
measurement strategies described in this work, namely, 
spatial postselection (Sec. Ill) and generahzed measure- 
ments (Sec. IV). Here, however, our discussion will be 
limited to the idealized case of point detectors. 

A. RSP of spatial qudits through spatial 

postselection 

The setup for spatial postselection of spatial qudits 
is exactly the same as that sketched in Fig. 2 with the 
double slit replaced by a D-slit array. A point detector 
constrained to move in the x direction between the focal 
and image planes of a lens will postselect the state [17] 

e 

\^{x,z)) = ipj{x,z)\j), (55) 
j=-e 

with the probability amplitude Lpj{x,z) given either by 
Eq. (4) in the Fresnel region or by Eq. (5) in the Fraun- 
hofer region. Thus if Alice wants to prepare Bob's qubit 
in the target state (54) she must find a position [x, z) to 
place her detector satisfying 

ipj{x, z) cx c*_y (56) 

It turns out that while for spatial qubits it is possi- 
ble to postselect arbitrary pure states through spatial 
postselection alone [17], the same is not true for spatial 
qudits [39]. This point, which is beyond the scope of the 
present work, will be addressed in more details elsewhere 
[40]. Therefore, when Alice shares the maximally entan- 
gled state (53) with Bob, she cannot prepare arbitrary 
pure states of his spatial qudit through this strategy. 

B. RSP of spatial qudits via generalized 

measurements 

To overcome the difficulties imposed by the spatial 
postselection for the RSP of spatial qudits, let us assume 
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that Alice implements the strategy sketched in Fig. 8 to 
measure her half of the maximally entangled pair (53). 
Again we assume that the overall two-photon state is, 
without loss of generality, |^)ab \H)a\H)b- Firstly, 
Alice implements the POVM shown in the dashed box of 
Fig. 8. Behind each slit there is a HWP followed by a 
phase shifter. These elements implement a unitary oper- 
ation which rotates the photon polarization by 26j and 
adds a phase shift (f>j, conditional to the passage of the 
photon through the slit j [30]. Note that this unitary 
operation differs from (24), for spatial qubits, only in the 
phase shift term. After coupling polarization and spatial 
degrees of freedom, the former is measured by a PBS. 
Secondly, Alice proceed with the measurement only if the 
outcome of the POVM was p = H; otherwise she does 
nothing. In case of success, the photon is detected by a 
single-pixel point detector in a fixed transverse position, 
a; = 0, at the focal plane of a lens, z = f. From Eqs. (5) 
and (55) one can check that this detection postselects the 
photon in a state proportional to J2j Therefore, after 
communicating the outcome of her measurement to Bob, 
his spatial qudit is left in the state 

e 

^ e'^^cos(20,)|-j). (57) 

In order to prepare the target state given by Eq. (54), 
she must set her HWPs and phase shifters, respectively, 
as follows: 

= ^cos-i(|c_j|), = arg(c_j)- (58) 

With the conditions (58), the outcome p = H oi the 
POVM implemented on Alice's photon occurs with prob- 
ability l/D for every setting of the plates. 

This strategy for RSP of spatial qudits differs from that 
for qubits (see Fig. 5) in three aspects: (i) the already 
mentioned unitary operation; (ii) just one output of the 
POVM is taken into account since for qudits the second 
output is useless [30]; (iii) it employs a single-pixel de- 
tector in a fixed transverse position at the focal plane. 
In this sense, this strategy is comparable to spatial post- 
selection in terms of the classical communication cost. 
However, RSP of spatial qudits via generalized measure- 
ment, as proposed here, offers major advantages from the 
practical point of view: it does not require the displace- 
ment of the detection system and may be experimentally 
implemented with programmable LCDs [27-29, 38]. This 
would result in a considerable reduction in experiment 
time, since instead of positioning the detector for each 
measurement, we would have only to update the prede- 
fined configuration (58) in the LCDs and measure the 
photon at the focal plane. More importantly, our strat- 
egy allows for the remote preparation of arbitrary pure 
states of spatial qudits which is not possible with spatial 
postselection alone. 



VII. CONCLUSION 



We have investigated, theoretically, the quantum com- 
munication protocol of RSP for pure states of spatial 
qubits. Two measurement strategies were studied, spa- 
tial postselection and a generalized measurement pro- 
posed by ourselves. In both cases we have analyzed the 
diffraction effects as well as the finite spatial resolution 
of the detection system over three figures of merit of the 
protocol, namely probability of preparation, fidelity, and 
purity of the remotely prepared state. It was shown that 
our strategy improves those three quantities compared 
with spatial postselection: it provides higher probabili- 
ties of preparation for the majority of the states in the 
Bloch sphere and allows their preparation with fideli- 
ties and purities close to the ideal value of one. Addi- 
tionally, we have presented a scheme for RSP of spatial 
qudits based on our measurement strategy proposed in 
[30], which allows for the preparation of arbitrary pure 
states unlike spatial postselection alone [39, 40]. In our 
schemes the crucial step in the realization of the proto- 
col, is the implementation of the POVM with the po- 
larization as the ancillary system. Experimentally, this 
could be achieved with programmable LCDs as recently 
demonstrated [28, 29, 38]. 

A possible extension of our work is the investigation 
of RSP of the angular spectrum of a single-photon field. 
This continuous degree of freedom allows one to encode 
a large amount of information. RSP could be a reliable 
method for transmitting this information. Recently, Wal- 
born et al. [41] proposed a QT scheme to accomplish such 
task for an unknown angular spectrum. In their scheme, 
Alice and Bob share a two-photon entangled state pro- 
duced by SPDC and a third photon, carrying the angu- 
lar spectrum to be transmitted, is in possession of Alice. 
Then, she performs a joint measurement on her photons 
by injecting them into a nonlinear up-conversion crystal 
cut for second harmonic generation and send the out- 
come to Bob through a classical channel. Depending on 
the outcome, Bob performs a unitary correction to re- 
cover the initial state. Nevertheless, there are two liming 
factors for the realization of the protocol. First, the state 
produced by SPDC is not equivalent to the maximally en- 
tangled Einstein-Podolsky-Rosen state [42] (which is re- 
quired for faithful teleportation) and thereby, the fidelity 
of the transmitted state is reduced. Second, the required 
nonlinear interaction in Alice's measurement strategy is 
very inefficient (^ 10^^), which compromises the overall 
efficiency of the scheme. If Alice had knowledge of the 
angular spectrum to be transmitted, RSP would be the 
better choice for this task, since the above difficulties im- 
posed by QT would be circumvented. First, because RSP 
could be faithfully performed with nonmaximally entan- 
gled states [6] and second, the measurement strategy em- 
ployed by Alice would not require nonlinear interactions. 
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